Vortex Core Excitations in Superconductors with Frustrated Antiferromagnetism 
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Motivated by recent discovery of cobalt oxide and organic superconductors, we apply an effective 
model with strong antiferromagnetic and superconducting pairing interaction to a related lattice 
structure. It is found that the antiferromagnetism is highly frustrated and a broken- time- reversal- 
symmetry chiral d + id'- wave pairing state prevails. In the mixed state, we have solved the local 
electronic structure near the vortex core and found no local induction of antiferromagnetism. This 
result is in striking contrast to the case of copper oxide superconductors. The calculated local 
density of states indicates the existence of low-lying quasiparticle bound states inside the vortex 
core, due to a fully gapped chiral pairing state. The prediction can be directly tested by scanning 
tunneling microscopy. 

PACS numbers: 74.70.-b, 74.20.Rp, 74.20.-z, 74.25.Jb 



Recently, superconductivity in the oxyhydrate 
Nao.35Co02 • I.3H2O below ~ 5K was discovered by 
Takata et al Eland confirmed immediately by several 
groups 0, 0, ETUI • Interesting features of this material 
include: (i) It is believed that the superconductivity 
comes from C0O2 layers, similar to that in copper-oxide 



cuprates. (ii) Co 4+ atoms have spin-| but form 
triangular lattice, which frustrates the antiferromag- 
netism (AF) and thus is a promising candidate for the 
occurrence of spin-liquid phases 6]. There are also 
organic superconductors like the ^-(BEDT-TTF^X 
materials which have a lattice structure very similar 
to the triangular lattice 0, 0. (in) Theoretically, the 
analysis based on the resonant valence bond the ory 
in the framework of the t-J model @, ES El 113 
or on the renormalizatio n g roup theory within the 
framework of t-U-J model indicates a wide window 
of broken-time-reversal-symmetry (BTRS) d+id'- wave 
pairing state in the phase diagram. Other theoretical 
groups 0, E3 EH proposed a p x + ipy-w&ve pairing state 
mediated by ferromagnetic fluctuations. Since the ferro- 
magnetism is insensitive to the detailed lattice structure, 
no frustration effect is expected on a triangular lattice. 
At this stage, we are not at a position to resolve the 
pairing state issue. Motivated by recent observation of 
a superconducting phase diagram of Na x Co02 • I.3H2O 
similar to that of the cuprate superconductors [jj, we 
consider in this Letter a spin-singlet pairing and study 
the nature of low-lying excitations around a vortex in 
these new superconductors with frustrated antiferro- 
magnetism. The results can be directly tested by further 
experiments such as scanning tunneling microscopy 
(STM), which likely will be carried out soon. 

In conventional s-wave superconductors such as NbSe2, 
the observed quasiparticle tunneling spectrum at the 
vortex core by Hess et al. E3 can be explained suc- 
cessfully in terms of the low-lying quasiparticle bound 
states as shown by Caroli, de Gennes, and Matricon [l^ . 
In copper-oxide cuprates, the condensate has a d-wave 















? 


r 








to- 











FIG. 1: Lattice structure of a superconductor with frus- 
trated antiferromagnetism. It has a nearest-neighbor hop- 
ping t and pairing interaction V on the bonds forming the 
two-dimensional lattice, a next-nearest-neighbor hopping in- 
tegral t' and V' along one diagonal of each plaquette, and an 
on-site Hubbard interaction U on each site. The sites A, B, 
and C labelled by filled circles are the vortex core center and 
its neighbors when the superconductor is in the mixed state. 



pairing symmetry. Theoretical study based on d-wave 
BCS model suggested 19] that, due to the existence of 
nodal quasiparticles, the local density of states (LDOS) 
at the d-wave vortex core exhibits a single broad peak 
at zero energy. However, the STM-measured local dif- 
ferential tunneling conductance at the vortex core cen- 
ter only exhibits a subgap double-peak structure 
or even no clear peak structure within the supercon- 
ducting gap [23|. The discrepancy between the ear- 
lier theoretical prediction and the experimental obser- 
vation stimulated various explanations. Recent inten- 
sive ex peri mental E3, E3, [2! EH E3 an d theoreti- 
cal [28 L l29l I30L l3lL 1321 l33j | studies seem to converge on 
an explanation in terms of the field-induced AF around 
the vortex core. When the AF is frustrated on a trian- 
gular lattice |3 EH ? one would expect a different nature 
of electronic excitations near the vortex. Previously, we 
have applied an effective microscopic mean-field model 
with competing AF and superconducting interactions to 
a square lattice, as relevant to the copper-oxide supercon- 
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ductors. This model generates very rich physics including 
the commensurate AF spin density wave ordering at un- 
doped systems, stripes at low doping, as well as the super- 
conducting states at optimal and overdoped regimes [2fj . 
Especially, within this model, it was found 29] that the 
AF ordering is induced around the vortex core, which 
explains several experimental observations on cuprates. 
Here we extend this effective model to a lattice as shown 
in Fig. m which interpolates between the square and tri- 
angular lattices. Analysis of this paper is also directly 
applicable to organic conductors that have such a lattice 
structure. 

The model consists of an on-site repulsion and off-site 
attraction. The former is solely responsible for the anti- 
ferromagnetism while the latter causes the superconduc- 
tivity. The mean- field Hamiltonian is written as: 

+ 5 (AyC *T C Jl + A ?i c ii c n) • (1) 

Here q ct annihilates an electron of spin <j at the ith site. 
The hopping integrals are respectively tij = t on the 
bonds forming the two-dimensional (2D) lattice and t' 
along one diagonal of each plaquette. The on-site repul- 
sion is U on each site. The quantities = (cj^c^), e^, 
and fi, are the electron density with spin a, the single site 
potential describing the scattering from impurities, and 
the chemical potential. The spin-singlet order parameter 
Aij = ^-{c^Cji — QjCjj) comes from the pairing inter- 
actions on the bond (V^- = V) and along one diagonal of 
each plaquette (Vij = V). The case of t' = and V = 
correspond to the model on a square lattice. With the 
application of an external magnetic field H, the Peierls 
phase factor is given by the integral ipij = ^ f* 1 A(r)-dr, 
where the superconducting flux quantum $0 = hc/2e and 
the vector potential A = x r in the symmetric gauge. 
The Hamiltonian can also be derived from the t-U-J 
model 0, |3tJ • We diagonalize Eq. JIJ by solving self- 
consistently the Bogoliubov-de Gennes equation: 
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subject to the self-consistency conditions for the electron 
density and the superconducting order parameter: = 
£„|< T | 2 /(£n) and n u = E„ |i#| 2 [l " /(£»)], and 

A« = ^E„K^r+ w a w "T) tanh (2tr)- Herethe 
quasiparticle wavefunction, corresponding to the eigen- 
value E n , consists of the component uf a for an elec- 
tron of spin g and the component vf a for a hole of op- 
posite spin <r. The single particle Hamiltonian reads 
Hij, a = —tije l< ^ ij + (Uriio- + ei — n)5ij. The Fermi distri- 
bution function is f(E) = l/[e E / kBT + 1]. Hereafter we 
measure the length in units of the lattice constant ao and 



the energy in units of t(> 0). As relevant to the supercon- 
ductivity in the cobalt oxides, we report results below for 
the case with t f = t and V = V. As a model calculation, 
we choose U = 4 and various values of V. We find that 
the antiferromagnetism is highly frustrated even when 
V = and the filling factor is one electron per site. Since 
the experiments on cobalt oxides were performed in the 
optimal or slightly overdoped regime 0, , we choose the 
filling factor rtf = G ni a /N x N y = 0.65, where N x , N y 
are the linear dimensions of the unit cell under consid- 
eration. The chosen band filling factor corresponds to 
an electron doping x = 0.35 We use an exact diag- 
onalization method to solve the BdG equation self- 
consistently. In zero field, the solution is found to be uni- 
form: There exits no AF spin density wave (SDW) and 
the superconducting bond order parameters exhibit the 
relation that = |A o |e"* , A y = |A |e^, A xy = |A | 
with = 27r/3, consistent with the results based on the 
t-J model (see, e.g., [l2^. We find typically |A | = 0.07 
and 0.15 for V = 2.5 and 3.0, respectively. This complex 
order parameter forms a broken-time-reversal-symmetry 
chiral pairing state. We then use this zero-field solution 
as the initial condition for iteration for the vortex prob- 
lem. 

When an external magnetic field is applied perpendic- 
ular to the 2D Co-Co plane, H = Hz (H cl <C H <C H c2 ), 
an Abrikosov vortex state is formed. Recent experiments 
suggest [sil EH that the cobalt-oxide superconductor is 
in the extreme type-II limit (the Ginzburg-Landau pa- 
rameter k > 100). Therefore, the screening effect from 
the supercurrent is negligible. By taking the strength of 
magnetic field, H = 2<$> /N x N y , such that the total flux 
enclosed is twice $0, the BdG equation Q is solved with 
the aid of the magnetic Bloch theorem |29j : 



v-k.,<j{T mn r) 



ik-R 



e »x(r,R)/2 UkCT ( r ) 

e -i X (r,R)/2 Vk)(T ( r ) 



(3) 



Here r is the position vector defined within a given unit 
cell, the vector R = mN x e x + nN y e y , k = e x + 

2ttI v 



wjk e v wltn l x,y 



0,1,...,M^- 



1 are the wavevectors 



defined in the first Brillouin zone of the vortex lattice, 
M X N X and M y N y are the linear dimension of the whole 
system, and the phase x( r > R) = |^ A(R) - r + 2mn7r. For 
the calculation, we consider a single magnetic unit cell of 
size N x x N y = 32 x 32. Typical results on the nature 
of the vortex core is displayed in Fig. [21 with V = 3. It 
shows that each unit cell accommodates two supercon- 
ducting vortices each carrying a flux quantum <£>o, which 
conforms to the above prescription for the magnetic field 
strength. The two vortices sit evenly on the diagonal 
of the unit cell, indicative of a diagonal-oriented square 
vortex. As shown in Fig.Efa)-(c), each component of the 
superconducting order parameter vanishes at the vortex 
core center and starts to increase to its bulk value, as one 
moves away from the core. Fig. Etd) displays the spa- 
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FIG. 3: The local density of states as a function of energy at 
three selected sites A (red line), B (green line), and C (blue 
line) for (a) V = 3 and (b) V = 2.5. The relative positions of 
the selected sites are indicated on Fig. [I] Also shown is the 
density of states for the zero- field uniform case. The energy 
is scaled to the maximum superconducting gap A macc . The 
other parameter values are the same as in Fig. except for 
T = 0.02. 



tial distribution of the staggered magnetization of the 
local SDW order as denned by Mi = (— l) 2 (n^ — n^). 
It is vanishingly small [A weak variation in the figure 
is due to the numerical accuracy]. The electron density 
n i — Xlcr n i° exhibits a Friedel oscillation around the 
vortex core (Fig. Efe)), in contrast to a monotonic de- 
crease of the electron density when a BCS d-wave vortex 
is approached. This is a direct manifestation of the vor- 
tex bound states inside the core. Finally, the vortex core 
shows a twofold symmetry, consistent with the underly- 
ing lattice structure under consideration. It is natural to 
expect that the vortex core on a triangular lattice would 
exhibit a sixfold symmetry. Most importantly, we find no 
local AF ordering around the vortex core in a supercon- 
ductor with a strong AF interaction, and therefore the 
vortex core is predicted to be a conventional one. Our 
results here are fundamentally different from the case in 
cuprates, where a local AF ordering is nucleated around 

the vortex core H3, E3, EE EE E3, ES El, EH . 
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FIG. 2: The three-dimensional display of the amplitude dis- 
tribution of (a-c) the three superconducting bond order pa- 
rameters, (d) the staggered magnetization M*, and (e) the 
electron density m = rii a in a magnetic unit cell contain- 
ing two vortices. The size of the unit cell is 32 x 32. Parameter 
values: *7 = 4, V = 3, n f = 0.65, and T = 0. 



The LDOS is defined by 

y k,n,cr 

+K"/f/'(£"' k + £)] , (4) 

where f'(E) is the derivative of the Fermi distribution 
function. The pi{E) is proportional to the local differ- 
ential tunneling conductance which could be measured 
by STM experiments. In Fig. [31 we plot the LDOS as a 
function of energy at the three selected sites around the 
vortex core center, which are labelled in Fig.Qfor various 
values of V. For comparison, we have also displayed the 
density of states for the zero-field uniform case, which 
is peaked at the maximum superconducting energy gap 
Amax- The coherent peaks at the gap edge are smeared 
out in the LDOS at the vortex core. Instead there are 
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FIG. 4: The spatial distribution of the LDOS at E = 
— 0.27A macc , at which the subgap resonance peak shows up 
in the LDOS at site C for V = 3.0. The other parameter 
values are the same as Fig.Hexcept for T = 0.02. 



netism. This investigation is relevant to recently discov- 
ered cobalt oxide and organic superconductors. We find 
no local induction of antiferromagnetism around the vor- 
tex core. Low-lying quasiparticle bound states are found 
inside the vortex core due to a fully gapped chiral pair- 
ing state. These core states can be directly measured 
by STM experiments regardless of the magnitude of the 
level spacing between them. To our knowledge, this is 
the first theoretical investigation of the nature of vortex 
core states in superconductors with frustrated antiferro- 
magnetism. 

We thank Y. Bang, M. Graf, Y. N. Joglekar, Z. Nussi- 
nov, and S. B. Shastry for useful discussions. This work 
was supported by the US Department of Energy. 



sharp subgap resonance peaks in the LDOS. Due to the 
energy level quantization of the bound states, the small 
gap around the Fermi energy follows approximately the 
relation, Ai ~ A^^/E^, similar to the s-wave case. 
For a strong pairing interaction, which leads to a large 
magnitude of the superconducting order parameter, the 
energy level spacing is wide (Fig.Efa)). For a weak pair- 
ing interaction, the level spacing is so small that only 
a thermally broadened peak is observed. This result is 
dramatically different from that in a BCS <i-wave super- 
conductor. In the latter, due to the existence of nodal 
quasiparticles, there always appears a broadened single 
peak around the Fermi energy regardless of how short 
the superconducting coherence length is. We propose to 
use the STM to investigate the core structure in these 
materials. For cobalt oxides, the superconducting energy 
would be very small, since T c ~ 5K at slightly over- 
doped regime. It is expected that a broadened single 
peak should be observed. It is a possibility that cobalt 
oxide superconductors have a spin-triplet pairing symme- 
try mediated by ferromagnetic fluctuations. One would 
expect therefore a competition between ferromagnetism 
and superconductivity in this case and vortex core exci- 
tations analogous to the ones in copper oxide cuprates 
may occur. This question is reserved for an interesting 
future study. 

We now look into the nature of the vortex states. In 
Fig. 01 we display the spatial distribution of the local 
density of states at a resonant energy level. It shows 
no long-range tails, in striking contrast to the case of 
d-wave superconductors, where tails runs along the gap 
nodal directions. Instead, the LDOS at the core level is 
trapped in an area characterized by the superconduct- 
ing coherence length, indicating that the core states are 
completely localized. 

In conclusion, we have studied the vortex core exci- 
tations in superconductors with frustrated antiferromag- 
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